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Abstract 



Quasiclassical approximation in the intrinsic description of the vortex filament dynamics is dis- 
cussed. Within this approximation the governing equations are given by elliptic system of quasi-linear 
PDEs of the first order. Dispersionless Da Rios system and dispersionless Hirota equation are among 
them. They describe motion of vortex filament with slow varying curvature and torsion without 
or with axial flow. Gradient catastrophe for governing equations is studied. It is shown that geo- 
metrically this catastrophe manifests as a fast oscillation of a filament curve around the rectifying 
plane which resembles the flutter of airfoils. Analytically it is the elliptic umbilic singularity in the 
terminology of the catastrophe theory. It is demonstrated that its double scaling regularization is 
QO . governed by the Painleve I equation. 
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1 Introduction 

Dynamics of vortices has been the subject of intensive study since the pioneering works of Helmholtz, 
Lord Kelvin and Poincare in the 19th century (see e.g. [UI2])- Relatively simple dynamics of a thin vortex 
filament in an incompressible inviscid fluid and in an infinite three-dimensional domain has attracted a 
particular interest. 

An approach based on the localized induction approximation (LIA) has been proposed in 1906 by 
Da Rios [3]. Within such approximation one assumes that the vorticity Q is confined to a thin vortex 
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core region and the core radius a of the vortex tube is much smaller than the radius of curvature R 
at the point and, hence, the induced velocity at this point is dominated by the contribution from the 
neighboring segment. Using the Biot-Savart formula 

= "J- / x -7 ^7 (!) 

4ttJ \x-x'\ z 

under the assumption that the circulation L of the vortex is a constant, one gets 

dX(s,t) r (L 



where X denotes a point on the filament, s is the arclenght, t is time, K is curvature, 6 denotes the 
binormal vector and L is the length of the segment of the filament whose contribution into the induced 
velocity is taken into account in the integral ([IJ . Assuming that one can neglect the variation of L along 
the filament, rescaling time by the constant In (— ) and using the expression of b in terms of X, one 
obtains 

9X - - 

— = Kb = X s x X ss , (3) 

Thus, within LIA approach, one has a simple binormal motion of the filament with the induced velocity 
proportional to the curvature K. Using the Serret-Frenet formulae for the space curves, one can transform 
(J3]) into the system of intrinsic equations for the curvature K and torsion r which is of the form [3] 

K t = -2K s r - Kt s , 

'K ss \ (4) 



r t = KK S - 2tt s 



K 



These LIA and Da Rios system have been rediscovered sixty years later in [31 [51 [6] . The detailed historical 
account of this story is presented in [7] . 

An important step in the development of the LIA approach has been done by Hasimoto in 1972 [8]. 
He observed that the system (0J is equivalent to the focusing nonlinear Schroedinger equation (NLS) 

#t + ^ ss + ^| 2 ^ = (5) 

via the transformation 

ip{s,t) =K{s,t)e i f ds ' T ( 8 '' t \ (6) 

The NLS equation was already widely known in plasma physics. Moreover in 1971 Zakharov and Shabat 
[9j discovered that NLS equation ([5]) is integrable by the inverse scattering transform (1ST) method. 
Hasimoto's result has demonstrated that the whole powerful machinery of the 1ST method (solitons, 
infinite sets of integrals of motion, symmetries etc.) is applicable to the vortex filament dynamics. This 
fact has lead to the explosion of interest to this field and since that time various aspects of the LIA in 
the vortex filament dynamics have been studied by different methods (see e.g. [TO]- [33]). Instabilities of 
vortex filament and formation of singularities were two important problems partially addressed during 
this period QS H3 [35] . 

At the same time it was clear that LIA is too restrictive since it neglects several important effects 
like axial flow, stretching, shortening, filament core structure etc. Attempts to go beyond LIA have been 
done in [37]- [42]. In particular, the inclusion of the axial flow into consideration has lead to the following 
expression for the induced velocity 



X t = ai K 2 t 



r 


H 


'2L\ 








+ a 2 


47T 







Kb + a 3 tx A + a 4 tx (Kb) s + Q (7) 



where t is the tangent vector, on are some parameters defined by the called matching condition, A is 
a certain vector and Q denotes the contribution from the nonlocal effects. In situation when one can 
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neglect Q and variation of — along the filament, the vector A can be chosen as (Kb) s . After a rescaling 
of time equation ([7]) takes the form [UJ 



= Kb + ll (^K 2 t + K s n + Krb 



(8) 



where n is the normal vector, 7! is a certain parameter related to the intensity of the axial flow. After 
the Hasimoto transformation equation ([5} becomes 

#t + V ss + + *7i + f MVs) = (9) 

that is the Hirota equation or the linear combination of the NLS and first higher NLS equation [43] . 
It is again an integrable equation. On the other hand it describes new effects such as axial flow, local 
stretching proportional to K s , binormal velocity depending on torsion r, effects which were not caught 
by LIA equation ©. It was shown in [5T] that axial flow affects essentially stability properties of the 
filament. The observation that the simple generalizations of equation §5§ or Da Rios system (£[]) may 
describe important nonlinear effects essential for vortex filament dynamics was one of the main results 
of [4TJ . Another observation made in [41] based on the paper [6] is that the study of the quasiclassical 
(or dispersionless) limit of the equation ([9]) clarifies essentially an analysis of destabilizing effects in the 
vortex filament motion. 

The latter result was strongly supported by the recent study of the modulation instability of the 
focusing NLS and higher NLS equations (see e.g. [H]-[SS]). In these papers it was shown that the 
dispersionless limit of the NLS equation (dNLS) is quite relevant to the study of the modulation instability. 
In this limit NLS is the elliptic system of quasilinear PDEs which exhibits the gradient catastrophe 
behavior at finite time. Behavior of solutions near to the points of catastrophe and beyond has been 
studied in details in [S5]-[S5]. It was argued in these papers that the gradient catastrophe for dNLS is a 
source of the modulation instability for NLS. 

The above two facts clearly indicate the importance of the integrable systems modeling certain effects 
of vortex filament dynamics and of their quasiclassical (or dispersionless) limit. 

Moreover, the results of the papers [3H1 HO], showing that the cutoff L/a is not a constant but 
strongly depends on curvature and torsion of filament, manifestly show the necessity to go beyond the 
approximation described by equations (|5]) and ([§]). Indeed, in general case, the parameters on in ([7]) are 
some functions of K, the vectors P and Q are certain superpositions of t, b, and n and then the induced 
velocity is of the form 

X t = A{K,r)t + B{K,T)n + C{K,T)b (10) 

where A,B,C are some functions. Equation (|10p implies certain intrinsic equations for K and r which, for 
particular choices of A,B,C, are integrable systems [TTJ different from ©. Relevance of these equations 
for the vortex filament dynamics is definitely of interest. 

Study of the quasiclassical approximation for the intrinsic equations associated with the motion (|10p 
and their implications for the vortex filament dynamics is the main goal of the present paper. Quasi- 
classical intrinsic systems describe motions of vortex filament with curvature and torsion slowly varying 
along the filament and in the motion, i.e. they are systems of equations for K — K(x,y) 1 t = r(x,y) 
where x and y are slow variables defined as x = es, y = et, e <C 1 (see e.g. [SS]). For the Da Rios system 
(HJ) one has 

K v = -2K x t-Kt x , 
T y = KK X - 2tt x 

or, in terms of Riemann invariants /3 = — r + iK 

Py = \{W+W*- (12) 
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For the intrinsic system associated with equation ([9]) (see [41]) and others the corresponding quasiclassical 
equations are of the form 

p y =X(p,p)P x (13) 



where the characteristic velocity is an appropriate function of P and /3. 
In this paper we will concentrate on the study of the equation 



|(3/3 + /3)+ 7i | (5/? 2 + m + /) +72 



P-/3 



Px (14) 



which is the the quasiclassical limit of the intrinsic equations corresponding to the evolution of the filament 
of the type (JTDJ) with 



A = 1i \k\ B = 11 K x+1 ^( 1 ^— 2 +KD\ , C = K + 1i Kt + 12 ( 1 JL- 2+ tD) ,1V, 



where D = K ln ^ T ^t^V^// ~~^ and 71,72 being real constants. First two terms in r.h.s. of ([T4"]) 
represents themselves the quasiclassical limit of equation ([9]) (see equations (3. 76), (3. 77)) from ([41] with 
71 = W), while the last term is of different type. Thought the corresponding contribution to the velocity 
v ([in]) is quite complicated (fT5|) the equation describing this flow separately is a very simple one. It is 
equivalent to the equation 

4> xx + (e% y = 0, cf = K 2 (16) 

which is the dispersionless Toda (dToda) or the Boyer-Finley equation [5n]-[B3]- dToda equation is 
integrable and plays an important role in various physical phenomena (see e.g. |63]-|67|). It is known 
also that dToda and dNLS equations belong to the same dToda/dNLS hierarchy (see e.g. [52] 167]). 
While the first two terms in r.h.s. of ([T4"]) represent "positive" members the last (dToda) term is the 
first "log- type" term. The presence of such dToda contribution into the induced velocity (fT0|) . even if 
not completely justified in a model independent way, could be of phenomenological importance since it 
might describe nonlinear effects different from those captured by equation ([9]). 

Ellipticity of equation (fT4|) or, in general, of (fl3]) suggests to write it, following to [68] [69] [70] , as the 
nonlinear Beltrami equation 

1 - i\{p,p) 

So, a solution of equation (fl"4"]l defines (see e.g. [55] [SSI HO] ) quasiconformal mapping if 

l + iX(P,P) 



l-i\(P,P) 



< 1. (18) 



In particular, any solution of the dDR system ([TT]) defines a quasiconformal mapping on the plane z. 
Quasiconformal mappings definitely play an important role in the vortex filament dynamics. 

According to the result of papers [50]- [58], equation ([14"]) well approximates the corresponding full 
equation of motion (equation at 72 = 0) until the derivatives P x and P y are not large. Analysis of 
behavior of solutions of equation ([M]) near the point of gradient catastrophe, where P x and P y become 
unbounded, is the main subject of our study. Some results concerning dDR equation (equation [T4l with 
7i = 72 = 0) have been reported briefly in [71) . Here, following the same approach, we will analyze the 
full equation (p~4|) . We will see that the inclusion of these two terms essentially modify the behavior of 
vortex filament. 

Our approach starts with the observation that the hodograph solutions of equation ([M]) describe 
critical points Wp = of the function W which obeys the Euler-Poisson-Darboux equation E (h, 5), i.e. 

(P-P)W^=l(We-Wj). (19) 

This fact essentially simplifies the analysis of the singular sector for equation ([T4"]) associated with the 
gradient catastrophe. 
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Since j3 x ,Py ~ ^ e P onl ^ s °f gradient catastrophe belong to a singular sector J\A sma in the space 
of solutions of the system (fT4|) for which Wpp — 0. Generic points of J\A smg are characterized by the 
condition Wppp ^ while for deeper singular subsectors J\A S ™ 9 , n > 2 all the derivatives of W vanish up 
to dp +1 W. For these subsectors one has 2n constraints on all variables parameterizing W. In particular, 
generic gradient catastrophe happens in a single point Xo,yo- 

Vortex filament exhibits a very particular behavior at the point of gradient catastrophe. While at a 
regular point (with bounded K, t,K x , t x ) a curve modeling the filament always lies of one side of the 
rectifying plane (spanned by the tangent vector t and the binormal vector b), at the point of gradient 
catastrophe it oscillates fast around the rectifying plane. In addition, center of the osculating sphere (and 
its radius) blows up to infinity. Moreover, the induced velocity V = X t becomes unbounded at the point 
Xq, 2/o- So, the gradient catastrophe for the system ([Tip give rises to the fast oscillations which resemble 
the classical airfoil flutter (see e.g. [72]). Numerical results [IS] show that such oscillations begin to 
expand along the filament |71] . Such a behavior of filament clearly indicates on its instability. 

Analytically at the point of gradient catastrophe one has a catastrophe of certain type. Using the 
double scaling limit technique we show that at the generic point of gradient catastrophe (i.e. in the sector 
Ai 1 ^ 719 ) one has the elliptic umbilic catastrophe in terminology of R. Thorn [73 with the principal part 
given by U 3 — 3UV 2 . Deeper gradient catastrophes are associated with higher singularities. For example, 
in the sector ng the corresponding singularity is given by Xg singularity from the V. Arnold list |74j 
with the principal part of the form U 4 — 6U 2 V 2 + V A . 

At the point of gradient catastrophe the quasiclassical approximation becomes invalid and the system 
(|14p should be substituted by its full dispersive version (by the intrinsic version of equation © at 
72 = 0). First order corrections can be obtained by the double scaling technique together with an 
appropriate modification of the function W* (j3*,(3 ) to a functional JW q (J3*,0 )dz*dz* such that the 
equation Wp = for critical points are substituted by the Euler-Lagrange equation = for W q . It 

is shown that for generic point of gradient catastrophe (i.e. the sector M.{ m9 ) the resulting equation for 
small corrections is equivalent to the Painleve I equation 

% = 6Q 2 - £. (20) 

The result of the paper [52] allows us to conjecture that any generic solution of the Da Rios system Q 
at the point of gradient catastrophe behaves as K = A"o + eK* and r = tq + et* where j3* = —t* + iK* 
is described by the tritronquee solution of the Painleve I equation. 

It is shown also that a class of double scaling regularizations of the higher gradient catastrophes 
associated with sectors A4^ ng is described by the equations 

n« = (n + 2)!n n+1 -f, n = 2,3,... . (21) 

The paper is organized as follows. In section 2 the relation of hodograph solutions with Euler-Poisson- 
Darboux equation is discussed. In section 3 it is shown that the hodograph solutions of dToda/dNLS 
equation describe critical points of a function. Hydrodynamic analogy and quasiconformal mappings for 
the governing equations are considered in section 4. Next section is devoted to an analysis of gradient 
catastrophe. Geometrical implications of gradient catastrophe and flutter of vortex filament are described 
in section 6. In section 7 it is shown that in the point of gradient catastrophe one has elliptic umbilic, 
Ag and other singularities. Double scaling regularization of these catastrophes is analyzed in section 8. 

2 Hodograph solutions and Euler-Poisson-Darboux equation 

The dDR system (|12p . equation (|14p and general two-component elliptic quasilinear system written 
in terms of Riemann invariants, i.e. (|13[) are linearizable by the classical hodograph transformation 
(x,y) -> (Re/3,Im/3) (see e.g. [75J[75]). One has 

x +Xy = O. (22) 

For equation (|13[) in the form of nonlinear Beltrami equation (1171) the corresponding linearized equation 
is 

(1 - i\)z p + (1 + iX)z/3 = 0. (23) 
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The compatibility condition for the system (|22j) is given by the second order equation 

(A - \)ypp = X-0p - Xpyj (24) 

solutions of which provide us with solutions of hodograph equations (|22j) . In the hyperbolic case, i.e. for 
hydrodynamic type equations, equation of the type (pMf has been widely used for the construction of the 
solutions and an analysis of their properties (see e.g. [7SJ [75] 177]). 

In our elliptic case, equation of type (j24j) written in term of variables v — —j3 — f3, u = — ((3 — /3) 2 /4 
and j3 defined as y — j3 v , i.e. the equation 

u(3 uu + 2f3 u + a(u)f3 uv = 0, a(u) > (25) 

has been applied to the study of the system 

u y + (uv) x = 0, Vy + vv x — a{u)u x = (26) 

which appears in the theory of wave processes in unstable media within nonlinear geometric optics 
approximation [78j [79]. At a{u) = 1 it is the quasiclassical limit of the focusing NLS equation (0 or 
the dDR system ([TTj) for u — K 2 , v = 2t. In particular, in the paper [SU], the behavior of the solutions 
for the system (|26[) near the singular point u = of equation ((25)) has been studied. For the hyperbolic 
version see also [HI] 

For the dDR system ([T2"]) equation (IM1) is particularly simple one 

(P-P)y P p = l(ye-yj)- (27) 

It is the Euler-Poisson-Darboux (EPD) equation of type E (|, |) (see e.g. [55]). For the dToda equation 
j3 y = j^pflx, equation ([24]) is of the form 



that is again the simple EPD equation of the type E (5,-5 

In contrast, for equation ([22]) with A = |(5/3 2 + 2(3/3 + (3 ) which corresponds to the second term in 
HH1) and to the quasiclassical limit of the first higher NLS equation (see formula ©), equation ([2~4"]) is 

(/3 2 - f)y m = 2(bp + 2(3)y p - 2(5/3 + 2p)yj. (29) 

For higher dNLS equations this equation is more and more complicated and hence an analysis similar to 
that performed in [791 81] becomes much more involved. 

A way to simplify drastically equations of type ([24l) associated with equation ([24]) can be extracted 
from the Tsarev's [83] discussion of the generalized hodograph equations. He noted that it is convenient 
to rewrite equation (|2"21 as 

{x + \y)j-y\j = 0. (30) 
Denoting x + \y = uj((3, (3), one observes that y = ^— = and equation (T5U|) assumes the form 

A — A 

oJs A- 

(31) 



lj — w A - A 

So, equation ([3H]) and, hence, the hodograph equation ([2"2l takes the form 

a; + AC8,^)i/-w09 J ^)=O (32) 

where u)(/3,p) is any function obeying the constraint (|31[) . In the form p2p the hodograph equations can 
be formulated also for multi-component hydrodynamic type systems [83] . 

An observation here is that for the dDR system (jTTj) and equation ([15]) one has 

Att 1 

(33) 



A- A 2{fi-P) 



(i 



and, hence, 

= 1 
w - 57 2(/3 - /3) ' 

Furthermore, these equations imply that 



(34) 



A = W>, u = Wp (35) 
where W"(/3, /3) and VF(/3, /3) are certain real valued functions. For instance, for equation 



|(3/3 2 + 2/3/3 + 3/3 2 ) + 7i^(5£ 3 + 3/3 2 /? + 3/3 2 /3 + 5/3 3 ) + 72 In f^r^ 



(36) 



In virtue of and ((SJ) the functions VF and W obey the EPD equation E (|, |) 



(37) 

; ( i r . 1 r- 1 

So, the hodograph equation (f3"0)l is of the form 

.x + yWjs - W = (38) 
where W and W are solutions of the same EPD equation E (h It remains to note that 



Wp-Wj w t 



Wp-Wj w m 



(39) 



is a solution of the equation (|24|) with A = Wp. 

Thus, in contrast to the usual approach in which hodograph solutions of the commuting quasilinear 
systems contributing into the r.h.s. of equation (fT4"|) can be found via different linear second order 
equations (see equations (fl"9)) . (j28j) . ([29])), in the approach described above all of them can be build via 
solutions of the same and simple EPD equation E (i, 5). 

EPD equations of different types are well studied and have number of remarkable properties (see e.g. 
[82"]). Their relevance for the theory of hydrodynamic type equations and Whitham equations has been 
understood in the papers |84j-|89j. 

3 Quasiclassical intrinsic equations for dToda/dNLS hierarchy 
and critical points of functions 

In addition to the observations made in the previous section, one readily concludes that the hodograph 
equation (|38|) can be represented as 

Wp = (40) 

where 

W = ~x(J3+Jfy+yW-W (41) 
and the function W obeys the EPD E (k, |) equation 

{P-?Wp~p = \(Wp-Wj). (42) 

Thus, one has 
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Proposition 3.1. Hodograph solutions of equations {1J$ describes critical points J^0| ) of the function W 
defined in (f^7| ) which obeys the EPD E (|, ^) equation. 

In the hyperbolic case, i.e. for the quasiclassical shallow water equation (or one-layer Benney system) 
this fact has been observed in [55J (see also |89|). 

So, in order to construct a solution of equation (TT4"|) one takes a function W of the form (|41l) with 

function W given by (|3l)j) . arbitrary functions W solving EPD A(|,i) equation and finds its critical 
points solving equation (|4T)|). The class of solutions obtained in this way is controlled by the class of 
solutions of the EPD E 5) equation. Various classes of solutions of this equation has been described 
a century ago (see e.g. 82J) as well as in recent papers concerning the application of the EPD equation 
in Whitham theory. For example, a class of solutions of the EPD E(e, e) equation is given by 



W = 



dX . x . / /3 



2ni 



(43) 



where S 1 is the positively oriented circle of large radius and u(X, t) :— X^fcLo ^ s an arbitrary polynomial 
of the complex variable A parameterized by the constants tk (see e.g. [55J 89 ). In particular, the function 
W (|4Tt with 72 = is given by (|43p with e = | and u = x + y(\ + 71 A 2 ). Parameters tk in u and W can 
be viewed as the variables parameterizing a class of initial data for equation (|14[) . 
In terms of curvature and torsion the EPD equation (|4"2"j) is of the form 



K(W KK + W TT ) + W K = 



(44) 



that is the axisymmetric three dimensional Laplace equation studied by Beltrami |90j (see also [75] ). 
Here A and r play the role of the radial and axial coordinates, respectively, in the cylindrical system of 
coordinates. 

Function W ED is 



W = —xt + y 
and hodograph equations are given by 

W K = y 



r- - \kA +7l {-t 3 + ^Kt 2 ) : - 2 hi K 



'2 K 



W K = 0, 



w 



(45) 



(46) 



W T = -x + y [2r + 71 (-3r 2 + 3Kt)] - W T = 0. 
The corresponding system of quasilinear PDEs in terms of K and r is of the form 





K x + 







K - 3rfrTK - ^ 
K 



K - 3 7 i - ^ 



A", 



2r - 7l ( Zt a - -if' 



r K =0, 
r K =0. 



(47) 



At 72 = it coincides with the dispersionless limit of the equations given in [41] . We note that the con- 
tribution from dToda terms (with 72 ^ 0) are quite different from those of dispersionless Hirota equation 
( ([47]) with 72 = 0) considered in [41]. Formally at small A, the last terms in (|47p becomes dominant. 
Essential difference between the dDR and dHirota type evolution and dToda contribution is manifested 
by the expression (TT5|) of the velocity of filament. The former are small for filaments close to a straight 
line (small A) and becomes larger for more curved filaments. In contrast, the dToda type contribution to 
the induced velocity (proportional to 72) is smaller for larger A and increase with decreasing A. Since in 
the approximation under consideration always A <C 1/a the dToda type contribution might be relevant 
for the description of some nonlinear and nonlocal effects appearing for slightly curved vortex filaments. 
In virtue of (14"2"]) at the critical points of the function W given by (|40[) or (|4"o| one has 



08 - fiWfiB = 



(48) 
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or 

K(W KK + W TT ) = 0. (49) 

The point K — is a singular point for these equations. Analysis of the behaviors of solutions of these 
equations near this point can be performed in a way similar to that discussed in 80J . 

Geometrically the case K = corresponds to a straight line. In this paper we will be interested in 
properties of curved vortex filament. In the rest of the paper we will assume that K ^ 0. So, equations 
(|48| and ([49]) take the form 

Wtf = (50) 

and 

W KK + W TT = 0. (51) 
We will see that these conditions simplify essentially an analysis of singular sector for equation (fI4"l) . 



4 Ellipticity of governing equations: hydrodynamic analogy and 
quasiconformal mappings 

Ellipticity of the dDR system (fTT|) and more general system (|47|) means that they have properties drasti- 
cally different from those of the hyperbolic version of the dDR system (or the system ([2"6"]> with a = — 1), 
that is the classical shallow water equation (see e.g. [76]). The latter has standard Ricmann wave solutions 
while in the elliptic case the Cauchy problem is ill posed. In spite of such a difference a representation 
of the elliptic systems in the form of hydrodynamic type systems, but with negative pressure, has been 
widely used in the analysis of their local properties (see [51 [751 179"! |4H |4"5] ) . 

Following this line one can rewrite the system (|4"T|) in terms of p = K 2 and v = 2 (l — j%A r — 
3 7 i (r 2 - \K 2 ) as 

Py + Ma = 0, 

p(Vy + VV X ) + P X = 

where 

P = ~\ (1 - 3 7 ir) 2 P 2 + ^7iV - 7172T (3p + 8r 3 ) + 2 72 (p + 2r 2 ) - 7 | U T - + hip) . (53) 

Here p, w, and P play the roles of density, velocity and pressure. Note that the system considered in 
[41] is not a system of quasi-linear equations. It contains the terms nonlinear in first order derivatives in 
addition to those quasi-linear ones in the system ([4"T]) with 72 = 0. Moreover our definition of the velocity 
and pressure is slightly different from those from [41] . Consequently, the system |52|) is composed exactly 
by the continuity and Euler equations in contrast to the system (3.79)-(3.80) from [4Tj . 

Properties of the system |52|) are fully determined by the properties of the pressure ([53]) . Within LIA 
the pressure P = — p 2 /2 is always negative and that causes the well-known instability of vortex filament. 
Inclusion of the Hirota and dToda contributions changes the situation drastically. For large curvature 
axial flow contributions dominate and the sign of the pressure is different for different values of K and r. 
For instance, if K > 2\/6|r| then the pressure P for large K and r is positive, while in the opposite case 
it is negative. Then, for the points with small r (near flattening points) one has stabilizing effect (P > 0) 
while, with large torsion r ~ a" 1 dHirota contribution (i.e. axial flow contribution) always destabilize 
the filament (recall that K ^ a -1 ). 

For small K (small density) pure dToda contribution equal to — 7! (in p + r 2 /p) dominates the pres- 
sure. So, for small K the sign of the pressure coincides with the sign of 4r 2 + 2K 2 \nK. Hence, once 
again one has a threshold dividing the regimes with positive and negative pressure and consequently 
stabilizing and destabilizing effects of contributions different from the original LIA model of Da Rios. If 
one considers the contribution only from axial flow (71 ^ 0, 72 = 0) the pressure is 

P = -\(1~^it) 2 P 2 + ^J 2 P 3 . (54) 
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P vanishes for K and r connected by the relation 37 2 A 6 = 8(1 — 3-fir) 2 K 4 , while it is positive for vortex 
filaments with curvature and torsion obeying the condition 



i-3r 

7i 



(55) 



and it is negative in the opposite case. 

This discussion shows that stability or instability of vortex filaments strongly depends on their geom- 
etry (relations between curvature and torsion). 

There is another way to formulate and emphasize an ellipticity of governing equations. The dDR 
system (|12[) and more generally system (|13l) are elliptic systems of quasilinear equations and, hence, 
under certain conditions their solutions define quasiconformal mappings on the plane (see e.g. [68 | 169 1 170]). 
Introducing the independent variable z — x + iy, one can rewrite general equation (|22|) in the form 



i + q(Ag) 



It is the nonlinear Beltrami equation [68, 69] 170) solution of which defines a quasiconformal mapping 

HA(/3,g) 



(z,C) — > iff the dilation coefficient /i = 1+ ^ x ^^ obeys the constraint \n\ < 1. This condition, i.e. 



H + *A(M| <1; 



|l-*A(/3,/3)| 

imposes constraints on /3 and, hence, on A" and r. For the dDR system ([TT]) . A = |(3/3 + /?) = — 2r + zA, 
and one has 

,o 4r 2 + (l-A) 2 , , 

Since K > one always has < 1 for any nonstraight filament. 

So, the slow motions of nonstraight vortex filaments within original LIA are described by quasicon- 
formal mappings on the plane defined by the equation 

2 + m + J) 

2-z(3/3 + /3y V ; 

For the generalization of the LIA given by equation ([9]) or (|14[) situation is more complicated. For equation 
(fl4| with 72 = one has 

, |2 8 + i [4(3/3 + g) + 3 7 i (5/3 2 + 2fl8 + g 2 )] 

Ml = = = =2- ( 60 ) 

8 - i [4(3/3 + /3) + 3 7 i(5/3 2 + 2/3/3 + /3 )] 
and the condition < 1| is verified if A and r and parameter 71 obey the constraint 

* ^ °' (2 + 2A) 2 + (4r - 6 7 ir) 2 + 3 7 iA (3 7 i A' - 8r) > °' (61) 

So, only part of the vortex filament configurations with axial flow is associated with the quasiconformal 
mappings, a sufficient condition being 71 t < 0. For the complete flow (|14|) the constraint on A, r, 71 and 
72 is more complicated. It is of interest to note that for pure dToda equation 

/3, = -0=? x (62) 



corresponding to the last term in (fl"4"]l . one has 
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So, solutions of the dToda equation define a quasiconformal mapping iff 72 < 0. At 72 = —1 equation 
(|62|) is equivalent to the equation 

^ + (e ) ro -O, = 21nX (64) 

and I /X I = and, hence, any solution of the dToda equation (IMl) defines a quasiconformal mapping 
on the plane. 

The results of the previous sections indicate also a relation between certain class of quasiconformal 
mappings on the plane and EPD equation. Indeed, solutions of equation (|4*tjf with A = Wp, i.e. solutions 
of the Beltrami equation 

*=±±5&A, (65) 
1-iWp 

are critical points Wp — of the function W = \x{(3 + /3) + yW — W which obeys the EPD E (|, i) 
equation. This means that the class of quasiconformal mappings on the plane defined by equation (|65p 
with arbitrary real- valued function W(/3,/3) obeying EPD E (4, -|) equation describe critical points of 
the function 

W = z(P + ~p~ -2iW) +z(/3 + /3 + 2iW) + W (66) 
where W is an arbitrary solution of the EPD E (h, |) equation. Inverse statement is valid too: critical 



points of the function (prJ)) obeying the EPD E (i |) equation and constraint 



1-iWja 



< 1 are described 



,2 ' 2; 

by quasiconformal mappings defined by equation (1661) . Similar results are valid also for hydrodynamic 
type system of e-type and other types. 

The interrelation between quasiconformal mappings and vortex filament dynamics will be studied 
elsewhere. 



5 Gradient catastrophe 

Hodograph equations Wp — and W-e — provide us with a single-valued solution of the dDR system 
(fT!?|) if the standard condition 

is satisfied. In virtue of the EPD equation (|42|) on the solution of the systems (jTTJ) and (fT2|) with K 7^ 
the function W obeys the equation = 0. Hence, at K 7^ 

A = \W PP \ 2 . (68) 

Consequently, the regular sector of equations (ITTI) and (|12l) for which hodograph equations are uniquely 
solvable is defined by the conditions 

Wp = 0, Wpp + 0. (69) 

In this sector both j3 and /3 x ,f3y are bounded. Indeed, differentiating the hodograph equation Wp = 
with respect to x and y and using the equation (|50)l . one gets 

* = -2^' *— ^ (7 °' 

Gradient catastrophe sector or singular sector for equations (JH]) is composed by those of their solutions 
for which j3 (i.e. curvature K and torsion r) remain bounded while derivatives j3 x , /3y explode (for standard 
definition see e.g. [76]). Thus, for the gradient catastrophe sector 

Wp = 0, Wpp - 0. (71) 

Condition of solvability of these equations is Wppp 7^ and, hence, generic singular sector is characterized 
by the conditions 

Wp = 0, Wpp=0, Wppp ^0. (72) 
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In terms of curvature and torsion these conditions at K ^ have the form 



W K =0, W r =0, W Kt = 0, Wkkk^O, W ttt ^0. (73) 



Equations Wkk = W TT = 0, Wkkt 0, Wktt are the consequences of (|73j) . 

Critical points of functions which obey the condition (|7T|) are called degenerate critical points (see e.g. 
[74] ) . Degenerate critical points are typical for families of functions parameterized by several variables 
[74] . In our case the family of functions W (|41l) is parameterized by x, y, 71, 72 and other variables defining 



initial values of (3 encoded in W. 

In the generic situation the last condition in ([?2"]) allow us to solve equations Wpp = 0, W-^g = with 

respect to (3 and j3, i.e. K and r. Substituting these expressions for K and r into the first equation (|72p. 
one gets two relations 

/i(^>y; 7i:72) = 0, M^, 2/571, 72) = (74) 

for equation (fl"4"]) . Thus, the generic gradient catastrophe for equation (| 14|) is characterized by two 
constraints on the independent variables x,y and parameters 71,72- For fixed 71 and 72 the gradient 

catastrophe for equation (|14[) occurs in a single point (xq, yo) for given W(f3, /3), i.e. for given initial data 
for K and t. For the dNLS equation 71 = 72 = this fact has been observed within another approach 
for the first time in [523 • In a different way it was derived in [SU] . 

For the quasiclassical counterpart of the Hirota equation ©, i.e. equation ([14")) with 72 = 0, one has a 
one-parametric family of points of gradient catastrophes for given initial data. Recall that the parameter 

71 in equation ([SJ and (|14[) encodes the matching data conditions for inner and outer expansions for the 
vortex filament with axial flow [37]- [41]. Thus, the meaning of a true physical "time" of the gradient 
catastrophe should be assigned to the minimal value 2/0(71) of the solutions 2/(71) of the equation ([74]) at 

72 = 0. It happens in the corresponding point #0(71) °f the filament. 

For equation (|14p and for the class of solutions of the dDR equation corresponding to a family of initial 

"data" W parameterized by additional variables, one can go deeper into the singular sector. Indeed it 
may happens that in addition to the conditions (|71j) in the singular sector the condition Wppp — is 
verified. For solutions for which Wpppp 7^ one has a subsector defined by 

Wp=0, Wpp=0, W m = 0, Wpppp^O. (75) 

Continuing this procedure one can show that the singular sector J\A sm9 of the dDR equation or equations 
of type (JUJ) has the following structure [551 [EH] 



M sing = [j M s n m9 (76) 

n>l 

where the subsector Aif^ 719 is defined as 

- d k W d n+2 W 

MT 9 = {(*, r, 7i; A P) G M smg : -j^r = 0, k = 1, 2, . . . , n + 1; £ 0}. (77) 

Solutions belonging to M.™ 19 are defined on a subspace of codimension 2n in the space of the variables 
x, y; 71, 72 given by equations 

/*(z,2/;7i,72)=0, fc = l,2,...,2n. (78) 

In order to have nonempty singular subsector M. sm9 one should have class of solutions characterized at 
least by 2n parameters. So, for equation ([T4"]) with fixed "initial data" W{j3, 0) one may have only two 
nonempty sectors M s ^ ng and Ai 6 ™ 3 ■ For appropriate family of the initial data one can have gradient 
catastrophe of arbitrary deepness. 

6 Flutter of filament 

A gradient catastrophe for the system ([47]) means a special behavior of a filament at the point xo at time 
yo. Using the formulae (fT0| and (fTB"]) . one observes that the induced velocity V has a particular behavior 
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at the point of gradient catastrophe. Within pure LIA 71 = 72 
bounded at the point xq. However, the acceleration given by 



0, V = Kb and, hence, velocity remain 



A = V t = 



(-2K s t - Kr s )b - KK s t — r 2 



(79) 



explodes at xq. 

Axial flow (71 7^ 0) and dToda contribution (72 7^ 0) qualitatively change this behavior. Indeed, the 
formula (fT5j) shows that at the point of gradient catastrophe the coefficient B, i.e. the normal component 
of the velocity becomes very large both for dHirota and dToda contributions. So, in the point of gradient 
catastrophe the local rate of self-stretching, proportional to X t ■ ft — B (see e.g. [HI [T7]) becomes very 
large. Note that such an effect is absent in pure dDR dynamics. Sign of normal component depends on 
the sign of K x and t x and, hence, at the point xq the filament begins to oscillate changing fast the sign 
of normal components of velocity. 

Local intrinsic consideration reveals a geometric origin and nature of such oscillations of the filament 
at the point of gradient catastrophe. It is well known, that the coordinates of a curve in a neighborhood 
of a point in the reference system formed by the tangent vector t, normal ft and binormal b and with 
origin in the point (see e.g. formula (6.7) from [91] and figure [lj are 



x 3 



K 2 3 
a - — s 3 +..., 
6 

— s 2 -s 3 + . . 

2 6 

Kt 



(80) 



1 

24 



(2K s t + Kt s )s 4 




X=X(s,t) 



Figure 1: Generic curve with r > at a regular point. 



At ordinary point where K, r and K x , t x are bounded the first terms in r.h.s. of (|80[) dominate and the 
curve, locally, is a twisted cubic. Its projections on the osculating plane (spanned by t and ft), normal 
plane (spanned by b and n) and rectifying plane (spanned by t and b) are parabola (12 — = 0), cusp 
(9Kxl - 2t 2 x\ = 0) and cubic (x 3 + = 0) respectively (fig. [2 see e.g. [91]). 

Characteristic feature of a curve (and a filament) in an ordinary point (i.e., outside the points of 
gradient catastrophe) is that it lies always on one side of positive direction of the normal and on one side 
of rectifying plane. 

At the point xq of gradient catastrophe the behavior of a curve changes drastically. Indeed, when K s 
and t s (i.e. eK x and er y ) become large, the r.h.s. of the formula (|80|) contains two parameters: small s 
and large K s and t s . So, the second terms in xi and X3 become relevant and different balance between 
s and K s (t s ) lead to different regimes in behavior of a curve around the origin. For instance, parabola 
in the osculating plane may change sign or even convert into cubic curve (fig. SJ while in the normal 
plane it could be a plane (3, 4) curve and so on (fig. [3]). So, around the point xq of gradient catastrophe 
a filament oscillates from one side of the rectifying plane to another and back. Such oscillation is quite 
similar to that of airfoil (see e.g. [72]) and one can refer to these oscillations of a filament in the point of 
gradient catastrophe as a flutter. 
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Figure 2: The figure show a curve with r > 0. At r < one should reflect the curve with respect to the 
axis t in the last figure. 




Figure 3: Different projections on the osculating plane of a vortex filament near to the gradient catas- 
trophe point. 



Another peculiarity of the point of gradient catastrophe can be seen in the behavior of the focal 
curve associated with filaments. The focal curve (see e.g. [HI]) is formed by the totality of centers of 
osculating spheres. Osculating sphere has a contact of third order with the curve in a points. Its center 
has coordinates 

X = X+±n-^b (81) 
where X is the position vector of a point on the curve. The radius R of the osculating sphere is given by 

2 K 2 T 2 +K 2 , N 

R2 = K 4 T 2 ' ■ (82) 

For a set of regular points the corresponding piece of the focal curve is regular too and it lies on a finite 
distance from the original curve. 

At the point of gradient catastrophe K s — > oo and, hence, the center of the osculating sphere as well 
as its radius blow up to infinity in the direction of the binormal vector (catastrophe direction, see figure 
[5]). So, the point of gradient catastrophe for the original curve at which both curvature and torsion are 
finite corresponds to a point of the focal curve at which the second focal curvature (—K s /K 2 t) becomes 
singular. We emphasize that even if our consideration for the original curve is purely local at the point 
of gradient catastrophe one observes a nonlocal effect of unbounded increase of the distance between a 
point on the original curve and the corresponding point on the associated focal curve. 

We note that the blow-up of Xo and R for the focal curve happens also at the flattening points at which 
t = (for a discussion on the flattening points see e.g. [92]). It would be of interest to understand fully 
a difference between the flattening points on a curve and points with unbounded K s in pure geometrical 
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Figure 4: Different projections on the normal plane of a vortex filament near to the gradient catastrophe 
point. 




terms. As an obvious difference one may observe that at the flattening point the projection of the curve 
on the osculating plane remains to be a parabola and the velocity V and acceleration V t remain bounded 
too. 



7 Elliptic umbilic and higher catastrophes 

In order to understand better the phenomenon described above one should analyze in more detail the 
structure of the singularity and the regularizing mechanism for the gradient catastrophe. For the focusing 
NLS equation such an analysis based on the e-expansion of the integrals of motions of the NLS/Toda 
equation has been performed in [52] . Here we will follow a different approach proposed recently in |93) . 

Thus, we consider a neighborhood of the gradient catastrophe point xo,yo and denote the values of 
j3 at this point by /3o- Following the double scaling limit method (see e.g. [94]- [95]) we will look for 
solutions of the system (|47p near the point of gradient catastrophe in the form 



X = Xq + € X , 

v = yo + e y , 
(3 = #> + 



(83) 



where e <C 1 and numbers a, er, 77 should be fixed by further consideration. In our approach the function 
W is the principal object for such an analysis. For the system (j47|) 



W = Wxx + W 2 y + W(P, /3) 



(84) 



where 



Wi 

w 2 



(3/? 2 + 20P + 3f) + ^ (5/3 3 + 3/3 2 £ + 3f(3 + hf) + 72 In (^A 
16 V 2i 



(85) 



For the solutions of the system (|47p one has 



Wp = w- 











(86) 



and also, for K 7^ 0, equation (|50|) , i.e. = 0. Furthermore differentiating EPD equation (|42l) with 

respect to (3, one has 

(87) 



2W p p + 2(J3-P)W pp p = W p 
Hence, at the point of gradient catastrophe fto, j3 due to equation (|50p and Wpp 









0. one gets 



W° - 

Y 0/3/3 



W°— 

0/3/3 



where f° :— f(Po,/3 ). So, at the expansion (|83)l for the function W, the third order derivatives W° 
and are the only nonvanishing third order derivatives in the singular sector M s , mg . 

0/30 jo 01 

As a result, an expansion of the function W near the point of gradient catastrophe is given by 



000 



W =W + e a W°x* + e a W°y* + \e a+ ^x*{fi* + /?*) + +e" +a ((W 2 °p + Wf 4 )P* + c.c) y* + 



(W 2 % + Wpp)f3* + {W 2 ° B - B + W^WP + C.C 



T : '>n 



(w^p/3* 



+ c.c. + 



(89) 



Consequently 



W . ^e a+ ^x* 



+ e^{W 2 p+W$)y* + 



F 3r/ 



2(W : ° 



100 



W$p)P* 



(W2°pp 



(90) 
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For generic point in M^ m9 the third term is smaller than the second one and the nontrivial balance of 
three terms in (|90f is achieved if a + ?/ = ?/ + a = 3ry, i.e. a = a = 2rj. 
Thus, near to the point of gradient catastrophe 

W = W +e 2r ' (W?x* + W$y*) + e 3ri W* (91) 

where 

W*(f3,(3*) = z*p* + z*~$~* + -a/3* 3 + -a~0* 3 (92) 

6 6 

and 

z* = \x* + (W 2 ° B + W$)y*, a = W$pp. (93) 
Denoting a 1 / 3 /?* = U + iV, one has 

W* = fU + gV + \u 3 - UV 2 (94) 

where / and g are certain linear combinations of x* and y*. In Thorn's catastrophe theory language a 
function of this form is known as describing elliptic umbilic catastrophe }73] . 

The type of behavior of f3 x and f3 y near to the point of gradient catastrophe follows from (|70p and 
Indeed, near to the point xg,yo, one has 

= Wp = e 3,1 W* B = e 2ri W* B , = z* + ^a(3* 2 . (95) 



= -£"" (~az*/2y 1/2 (96) 



Hence f3* = (-^f^' 2 and 

P» 

The elliptic umbilic singularity (|94|) and formula f|96f) describe generic point of gradient catastrophe 
which corresponds to the sector AA s ^ ng . For variable parameters 71, 72 and other parameters describing 
the family of initial data, higher singularities corresponding to the singular sectors A4^ ng discussed in 
section 5 may appear. 

For instance, if for some 71, 72, Wggg = and Wq SBB ^ then instead of the expansions one 

has 

W =W + e a W^x* + e a W%y* + ^e a+v x*(/3* + /f ) + e 7 ' +a ([W 2 ° B + W°)f3* + c.c.) y* + 
+ — [(W2°pp + W° B3 )P* 2 + {W 2 y g + + c.c] y* + — [WjipppP 



4 + c.c. I + 



So, 



Wp* =\e a *>x* + s^(W 2 ° B + W° B )y* + ^ [2(W 2 °pp + W° BB W + (W 2 ^ + Wfetf 
+ -^W°pppp0* 3 + --- = O. 



(97) 



(98) 



and, hence, in this case 

a = a = 377. (99) 



W = W + e 3n (W°x* + W°y*) + e ir 'W* + ... (100) 

where 

W* = z*(3* + z*f + —a/3* 4 + —hp* 4 (101) 
24 24 
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and z* = \x* + (W 2 p + W%)y* as in the previous case (when ^= 0), while now a — Wpppp. 

Denoting 5 1 / 4 /?* = U + iV, one gets 

W* ^lU + mV+^U i -^U 2 V 2 + ^V 4 . (102) 

The last three terms reproduce the canonical form of a (nonsimple) unimodular parabolic singularity Xg 
in terms of Arnold's classification. 

Near to the point of gradient catastrophe xo , yo the behavior of /3* is given by 

Q = Wp = Wp* =z*+ (103) 



Therefore /3* = f-^-J and 



If 6 n 1/3 



«- = 3l-^l • (1 ° 4) 



In the sector Ai^" 19 generically one has a = a = (n + 2)rj and therefore 

W = W + £ ( n+2) "W* (105) 

where 

(n + 2)! (n + 2)! 

with a„ = (d r > +2 W)° and 

~ z *-"/(n+i)_ (107) 

Appearance of these singularities and particularly the time order of their appearance strongly depends 
on concrete structure and properties of vortex filament. 

For nongeneric point of singular sector one might have a different behavior of function W. For instance, 
for the singular sector M 6 ™ 3 in situation when 

Wipifa ,£„) + fWoJo) - 0, (108) 

the third term in (|89p of the order 2?7 + a becomes relevant and the balance of order requires that 
a + T) = a + 2?y = 3f?. Hence, in such situation a — 2rj, a — rj and instead of the expression (|90[) one has 

W =W + e^W^y* + £ 2 "i(/3 +P )x* + e^W* . (109) 



where 

W* =\x*{P* + \ [(W 2 % + W° fj )(3* 2 + {W 2 % + W%)P*T + cc. 



(110) 



So, in the case (|108[) the variables x* and y* play a quite different role in contrast to the generic case 
(ji?2"j) . Such type of behavior will be studied elsewhere. 



8 Double scaling regularization. Painleve I and other equations 

At the point of gradient catastrophe the quasiclassical approximation ceases to be valid and contribution 
of the higher order derivatives terms in the governing equation becomes relevant. A standard way to 
deal with this problem is to begin with the full (dispersionfull) version of the governing equation (i.e. Da 
Rios system (j4]) or intrinsic version of equation (j9))) and, then, performing an appropriate double scaling 
limit, calculate the required regularization terms. For the NLS equation ([5]) such an analysis has been 
performed in |52) . 
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Here we will follow another approach discussed in [931 171] . The basic idea is to incorporate higher 
order derivatives terms without any reference to the exact (dispersionfull) equation. 

The first step is based on the observation (section 3) that governing equations describe critical points 
of the function W. Thus, in order to incorporate derivatives of p, it is quite natural to modify W to 
W q (/3, P, P x , P x , (3 y , f3 y ) in such a way that equation (|40|) is substituted by the Euler-Lagrange equations 



dWg 

dp 



= 0, 



^=0 
dp 



(111) 



for the functional J W q dzdz. 

Second point is to use the specific properties of the function W . Namely, the real-valuedness of W 



W{p,p) = W{p,f3) 



(112) 



and the fact that in the perturbations W* (see formulas (j9"2"]) . (I106p ) the contributions of j3* and /3* are 
separated. It is natural to look for the modifications of W which preserve these properties. Thus, near 
the point of gradient catastrophe one has the following W q (omitting inessential terms in (|9"Tj) and (|92j) 
and choosing rj = 1) 



W q = Wo 



bey 



* 2 

-6)8- 



(113) 



where S and b are appropriate constants. Euler-Lagrange equations 

-0, 

= 



dW q , 


'dW q 


d/3* 1 


ydp* z , 


dW q 


( dW q 


— * 

dp 





(114) 



implies that 6 = 3. So, one has the equation 

1 

2 



a/3* 



0. 



This equation is converted into the classical Painleve I equation 

n a = 6n 2 - c 

by the change of variables 

z* = \i, 



^ = 4b n 



(115) 

(116) 
(117) 



with a A 5 = -246 2 . 

So, the regularized behavior at the point of the gradient catastrophe for the vortex filament dynamics 
is governed by the Painleve I equation. The relevance of this equation for the NLS/Toda system near 
to the critical point has been first observed in a different way in [52] . The correspondence between the 
solution of equation (|115[) and those given by the formula (5.21) in [52] is (at 6 = 1) u = K 2 , v = 2t and 



4/5 



(re**) 



-1/5 



36K 



1/5 



f3 = - ibe 2/5 (r 



36K 



3/5 



(118) 



where x, P, r and tp are defined in (2.19) and (2.22) of [55]. One also concludes that e = e 2 / 5 . 

In the paper ([55]) it was conjectured that any generic solution of the NLS equation near to the critical 
point is described by the tritronquee solution of the Painleve I equation. Extension of this observation 
to our case allows us to formulate the 
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Conjecture 8.1. Generic behavior of the vortex filament near the point 3:0,2/0 of gradient catastrophe 
within the localized induction approximation is given by 

*.,«) « a - • (-H) a, 1/5 ^) (us) 

where fio *s iAe tritronquee solution of the Painleve I equation ill 6}) . 

Such an approach can be extended also to the singular sectors M.^ ng . In these sectors the modified 
W q , for example, can be chosen in the following form 

w qn = Wo + e n (^( z *[3* + rT) + ^2y a ^* n+2 + (^T2)i^* ,i+2 ) +£Sn (\ bn ^ + \ v ^* 2 *> 

(120) 

where 6 n are appropriate constants. The Euler-Lagrange equation readily implies that 5 n — n. So, the 
modified W is given by 

Wqn = W + e n fyz*[3* +??) + ^-2_ a „/r" +2 + ^J-^a^* n+2 + \b n [3* 2 z , + \kTIA (121) 
and the corresponding Euler-Lagrange equation for ft* is 

5^151^' + f (122) 

By the appropriate change of variables 

z* = A£, j 8* = _A_n, 2 a „A 3n+2 = -(-26„) n+1 (n+l)!(n + 2)! (123) 
one gets the equation 

% = (n + 2)!fr +1 -f. (124) 

Eventual relation of these equations with appropriate limits of NLS hierarchy, higher Painleve equations 
and other types of regularizations containing higher order derivatives will be considered in a separate 
paper. 
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